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We present model independent predictions on the short range behavior of the energies of the gluonic excitations 
between static quarks (hybrid potentials) in an effective field theory framework (pNRQCD). 
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1. INTRODUCTION 

One of the clearest way of studying the heavy 
hybrid systems is from the energies of the gluonic 
excitations between static quarks, the so called 
excited gluonic potentials Vh ■ These energies are 
calculated on the lattice using Wilson loops with 
insertions of operators of a definite symmetry (see 
below) in the end point strings. The frame is 
supplied by the Born-Oppenheimer approxima- 
tion. In this picture the quarks and the gluons 
are respectively slow and fast degrees of freedom 
and in the limit Aqcd/*™ — * oo, (m = mass of 
the quark), the static gluonic energies Vh are 
treated as adiabatic potentials for the QQ pair 
and inserted in the Schrodinger equation to get 
the masses of the hybrids. Of course, the shape 
of Vh as a function of the QQ distance r is an 
interesting observable of QCD. 

Here, we show how the appropriate QCD effec- 
tive field theory description of the heavy quark 
bound systems can provide us with model inde- 
pendent predictions on the short range behavior 
of Vff. 

2. pNRQCD 

The existence of a hierarchy of energy scales 
m ^ mv 3> mv^ , v <^ 1 being the quark veloc- 
ity, in the heavy quark non-relativistic bound sys- 
tems allows us to substitute scale by scale QCD 
with simpler but equivalent effective theories. In- 
tegrating out the hard scale m one passes from 
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QCD to NonRelativisticQCD (NRQCD), while 
integrating out the soft scale mv one passes from 
NRQCD to potential NRQCD (pNRQCD) §]. In 
this last effective theory only ultrasoft (US=with 
energy much smaller than mv) degrees of freedom 
remain dynamical. The surviving fields are the 
Q-Q states (with US energy) and the US gluons. 
The Q-Q states can be decomposed into a singlet 
(S) and an octet (O) under color transformation. 
The relative coordinate r = Xi — X2, whose typi- 
cal size is the inverse of the soft scale, is explicit 
and can be considered as small with respect to 
the remaining (US) dynamical lengths in the sys- 
tem. Hence the gluon fields can be systematically 
expanded in r (multipole expansion). Therefore 
the pNRQCD Lagrangian is constructed not only 
order by order in 1/m, but also order by order in 
r. As a typical feature of an effective theory, all 
the non-analytic behavior in r is encoded in the 
matching coefficients, which can be interpreted as 
potential-like terms. The pNRQCD Lagrangian 
density is given at the next to leading order in 
the multipole expansion by: 
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where R = (xi +X2)/2, S = S(r,R,<) and O = 
0(r, R, t) are the singlet and octet wave functions 
respectively. All the gauge fields in Eq. (h^) are 



evaluated in R and t. In particular E = E(R, t) 
and iDoO = idoO - glAoCR, t),0]. Vs and Vo are 
the singlet and octet matching potentials; Va and 
Vb are matching coefficients. 

In QCD another physically relevant scale has to 
be considered, i.e Aqcd, the scale at which the 
nonperturbative effects become relevant. As far 
as mv :s> Aqcd, the matching between NRQCD 
and pNRQCD can be performed perturbatively. 

We believe that pNRQCD is the most appro- 
priate effective theory to describe the physics of 
heavy quark bound systems. For details see 

3. GLUELUMPS in pNRQCD [g] 

Although pNRQCD is originally designed to 
study Q-Q systems of large but finite mass, it is 
interesting to study its static limit. We shall re- 
strict ourselves to the case of pure gluodynamics. 
We call gluelump the adjoint source in presence 
of a gluonic field 

H{-R,r,t) ^ H''{R,t)0''(R,r,t). (2) 

In the static limit there are no space derivatives 
in the Lagrangian, and hence r and R are good 
quantum numbers. The spectrum then consists of 
static energies which depend on r (translation in- 
variance forbids R dependences) and on the only 
other scale in this problem, Aqcd- 

In the limit Aqcd ^1/?", the spectrum of the 
theory can be read from the Lagrangian (|l|). In 
particular, the leading order solution corresponds 
to the zeroth order of the multipole expansion. At 
this order the dynamics of the singlet and octet 
fields decouple. Hence, the gluonic excitations 
between static quarks in the short-distance limit 
correspond to the gluelumps. Depending on the 
glue operator H and its symmetries, the gluelump 
operator (H) describes a specific gluonic excitation 
between static quarks and its static energy Vh- 

In NRQCD (as in QCD) the gluonic excitations 
between static quarks have the same symmetries 
of a diatomic molecule plus charge conjugation. 
In the centre-of-mass system these correspond to 
the symmetry group Dooh (substituting the par- 
ity generator by CP). According to it the mass 
eigenstates are classified in terms of the angu- 
lar momentum along the quark-antiquark axes 



{\Lz\ — 0,1,2,... to which one gives the tradi- 
tional names E,n, A, ...), CP (even, g, or odd, 
w), and the refiection properties with respect to a 
plane passing through the quark-antiquark axes 
(even, +, or odd, — ). Only the S states are not 
degenerate with respect to the reflection symme- 
try. In Fig. 1 we report the lattice measurement 
of Vh obtained considering Wilson loops with op- 
erators of the appropriate symmetry inserted at 
the end points 

In pNRQCD at lowest order in the multipole 
expansion, besides the already mentioned symme- 
tries, extra symmetries for the gluonic excitation 
between static quarks appear. The glue dynam- 
ics no longer involves the relative coordinate r. 
Therefore, the glue associated with a gluonic exci- 
tation between static quarks acquires a spherical 
symmetry. In the centre-of-mass system gluonic 
excitations between static quarks are, therefore, 
classified according to representations of 0(3)x 
C, which we summarize by L, the angular mo- 
mentum, CP and refiection with respect a plane 
passing through the quark-antiquark axes. Since 
this group is larger than that one of NRQCD, 
several gluonic excitations between static quarks 
are expected to be approximately degenerate in 
pNRQCD, i.e. in the short-distance limit r <^ 
1/Aqcd- We illustrate this point by building up 
all operators, H, up to dimension 3 and by clas- 
sifying them according to their quantum num- 
bers in NRQCD and pNRQCD in Tab. In 
Tab. all the operators are intended evaluated 
in the centre-of-mass coordinates. is not dis- 
played since it corresponds to the singlet state. 
The prime indicates excited states of the same 
quantum numbers. The chosen operators for the 
n and A states are not eigenstates of the refiec- 
tion operator. This is not important since these 
states are degenerate with respect to this symme- 
try. The operators E and D x B (B and D x E) 
have the same quantum numbers and are related 
by the equations of motion so that they project 
over the same states. From the results of Tab. 
|l| the following degeneracies are expected in the 
short-distance limit: 

s+'-n, ; s;~n;^A,; 

^ n„ ; j:+^k^a^. (3) 



Table 1 

Matching operators H for the E, 11 and A gluonic excitations between static quarks in pNRQCD up to 
dimensions 3. The covariant derivative is understood in the adjoint representation. D • B and D • E do 
not appear, the first because it is identically zero after using the Jacobi identity, while the second gives 
vanishing contributions after using the equations of motion. 





Gluelumps 
(Tr{Oaffn) 


L = 1 




L = 2 




r • E ,r ■ (D X B) 












(r.D)(r.B) 




r X E ,r X (D X B) 






n; 






r X ((r •D)B + D(r -B)) 


A, 




(r X D)*(r 


X By + (r X Dp(r x B)' 


S+ 






(r.D)(r.E) 




r-B ,r • (D X E) 








r X B ,r X (D X E) 






K 






r X ((r •D)E + D(r-E)) 






(r X D)*(r 


X Ep + (r X Dp(r x E)' 



Similar observations have also been made in 
In pNRQCD they emerge in a quite clear and 
straightforward way. Moreover, here we can write 
explicitly the relevant operators. 

So far we have just used the symmetries of pN- 
RQCD at lowest order in the multipole expansion. 
In fact we can go beyond that and predict the 
shape of the static energies by actually studying 
the correlators 

(0liI(R,r,T/2)i?t(R',r',-T/2)|0} ~ 

,53(R-R')^'(r-r')e-'^^"M (4) 

for large T. At leading order in the multipole 
expansion we obtain 



VH{r) = Vo{t) + 



(5) 



+ - ln(i/"(T/2)0(r/2, ~T/2T:^H\~TI2))), 

where the T — > oo limit is understood. The gen- 
eral structure of the gluonic correlator is the fol- 
lowing (the contribution from the continuum is 
included in the dots) 



(iJ'^(T/2)0(r/2, -T/2)lfH''{^T/2)) 



\ nonpert. 
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Since we are in the static limit, l/T <^ Aqcd 
Ah < A^ < . . ., one can approximate the right- 
hand side of Eq. (||) for T — > oo by just keeping 
the first exponential of Eq. (^. Then we get at 
leading order in the multipole expansion 

VH{r) = Vo{r)+AH. (7) 

Formula states that at leading order in the 
multipole expansion the short-distance behavior 
of the static energies for the gluonic excitations 
between static quarks is described by the pertur- 
bative octet potential plus a nonperturbative con- 
stant. The constant Ah depends in general on 
the particular operator H, i.e. on the particular 
gluonic excitation between static quarks. Ah is 
the same for operators identifying states which 
are degenerate. Notice also that Eq. (0) can be 
systematically improved by calculating higher or- 
ders in the multipole expansion. In particular, 
one can look at how the 0(3)xC symmetry is 
softly broken to Dooh in the short-distance limit. 
Let us compare our results against the best avail- 
able lattice data (see Fig. |l| and First of all 
we observe a tendency of all the static energies ex- 
cept Hg to go up at short-distances which may be 
interpreted as an indication that they want to fol- 
low the octet potential shape as required by (^. 




Figure 1. Energies Vh for different gluonic excita- 
tion between static quarks at distance r from the 
quenched lattice measurements of Q], rg ~ 0.5 
fm. The picture is taken from [Q. 

The Ilu — ^'^d Au — Tl'^ — Si" static energies 
also show a strong tendency to form a degenerate 
doublet and triplet respectively at short-distances 
as precisely predicted in (|^). However, something 
strange is observed for the g static energies. Ac- 
cording to (^) a degenerate doublet and triplet 
should be observed. We can clearly see the dou- 
blet Ug — E+' but, as mentioned before, the shape 
of Ug is not compatible with the octet potential. 
On the other hand Ag — E~ miss a state 11^ to 
complete the triplet. We suspect that the plot- 
ted Ug is a superposition of the real Ug in the 
doublet and the in the triplet. This would 
explain its peculiar behavior at small r. We ex- 
pect that future lattice simulations will be able 
to disentangle these two states and confirm our 
predictions. It is interesting to notice that also 
the hierarchy of the states, as displayed in Fig. ||, 
is reflected in the dimensionality of the operators 
of Table |^. Lower lying states are characterized 
by lower dimensional operators if we use a min- 
imal bases with no electric fields. In particular 
this allows to understand on simple dimensional 
grounds the highly non-trivial fact that the dou- 
blet and triplet static energies lay close together 
for the g states but quite far apart for the u ones. 

The results given above allow us to relate, in 
the short-distance limit, the behavior of the en- 



ergies for the gluonic excitation between static 
quarks with the large time behavior of some 
gluonic correlators, in particular with their cor- 
relation length. It is particularly appealing 
that we can extract results for the gauge in- 
variant two-point correlator for the gluon field 
strength tensor: {0\F^Jtmt,0)lfF'^,{0)\0),^ 
being the adjoint string. One can pa- 
rameterize this correlator as a function of 
two scalar functions: {0\E,''{t)(l){t,0)l'l^'E\0)\0) 
and {0\B''{t)(l)(t,0)l'l^B''{0)\0) with correlations 
lengths: Te — 1/Ae and Tb = l/A^ respec- 
tively. From the results of Ref. displayed in 
Fig. H we can conclude that 

Te<Tb. (8) 

So far lattice simulations of the gauge invariant 
two-point correlator for the gluon field strength 
tensor have not reached enough precision to con- 
firm this behavior . This would be a nice cross- 
check of lattice simulations. On the other hand, 
recently, a sum rule calculation found evidence 
in favor of Eq. and also a lattice computation 
of the gluclump masses ^ seem to confirm Eq. 
(^. It would be highly desirable to have more 
precise lattice data in the short-distance limit in 
order to test our results more quantitatively. 
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